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ABSTRACT 

The intrinsic alignment (IA) of galaxies has been shown to be a significant barrier to 
. precision cosmic shear measurements. (Zhang, 2010, ApJ, 720, 1090) proposed a self- 

04 ] calibration technique for the power spectrum to calculate the induced gravitational 

shear-galaxy intrinsic ellipticity correlation (GI) in weak lensing surveys with photo-z 
measurements which is expected to reduce the IA contamination by at least a factor of 
10 for currently proposed surveys. We confirm this using an independent analysis and 
propose an expansion to the self-calibration technique for the bispectrum in order to 
• calculate the dominant IA gravitational shear-gravitational shear-intrinsic ellipticity 

correlation (GGI) contamination. We first establish an estimator to extract the galaxy 
Q . density-density-intrinsic ellipticity (ggl) correlation from the galaxy ellipticity-density- 

density measurement for a photo-z galaxy sample. We then develop a relation between 
the GGI and ggl bispectra, which allows for the estimation and removal of the GGI 
correlation from the cosmic shear signal. We explore the performance of these two 
methods, compare to other possible sources of error, and show that the GGI self- 
calibration technique can potentially reduce the IA contamination by up to a factor of 
5-10 for all but a few bin choices, thus reducing the contamination to the percent level. 
The self-calibration is less accurate for adjacent bins, but still allows for a factor of 
QO ■ three reduction in the IA contamination. The self-calibration thus promises to be an 

efficient technique to isolate both the 2-point and 3-point intrinsic alignment signals 
. from weak lensing measurements. 

Key words: gravitational lensing - cosmology 



1 INTRODUCTION 

Weak gravitational lensing due to large scale structure (cosmic shear) has emerged as a powerful cosmological probe in order 
to map the distribution of dark matter in the universe and to characterise the equation of state of dark energy, improving 
constraints on the equation of state of dark energy and the matter fluctuation amplitude parameter by factors of 2 to 4 (see 
for example Bacon, Refregier & Ellis (2000); Brown et al. (2003); Eisenstein, Hu & Tegmark (1999); Fu, Wu & Yu (2009); 
Hoekstra et al. (2002); Hu & Tegmark (1999); Hu (2002); Jarvis et al. (2003); Joudaki, Cooray & Holz (2009); Massey 
et al. (2005); Pen et al. (2003); Rhodes, Refregier & Groth (2001); Schrabback et al. (2010); Van Waerbeke et al. (2000, 
2002); Zaldarriaga, Spergel & Seljak (1997) and references therein.) Gravitational lensing has also been shown to be very 
useful to test the nature of gravity at cosmological distance scales (see for example the partial list Acquaviva et al. (2008); 
Bean & Tangmatitham (2010); Capozziello, Cardone & Troisi (2006); Daniel et al. (2008, 2010); Dossett, Moldenhauer & 
Ishak (2011); Huterer & Linder (2007); Ishak, Upadhye & Spergel (2006); Ishak & Dossett (2009); Linder & Cahn (2007); 
Schmidt (2008); Song (2005); Thomas, Abdalla & Weller (2009); Toreno, Semboloni & Schrabback (2010); Zhang et al. 
(2007); Zhao et al. (2006, 2009).) 

In addition to the constraints obtained from the 2-point cosmic shear correlation and the corresponding shear power 
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spectrum, the 3-point cosmic shear correlation and the shear bispectrum have been shown to break degeneracies in the 
cosmological parameters that the power spectrum alone does not (Takada & Jain 2003; Vafaei et al. 2010). For example, 
the results of Takada & Jain (2004) showed that a deep lensing survey should be able to improve the constraints on the dark 
energy parameters and the matter fluctuation amplitude by a further factor of 2-3 using the bispectrum, and most recently, 
Semboloni et al. (2010) derived parameter constraints by measuring the third order moment of the aperture mass measure 
using weak lensing data from the HST COSMOS survey. They found independent results consistent with WMAP7 best-fit 
cosmology, and an improved constraint when combined with the 2-point correlation. Ongoing, future and proposed lensing 
surveys (e.g. CFHTLS 1 , DES 2 , EUCLID 3 , HSC 4 , HST 5 , JWST 6 , LSST 7 , Pan-STARRS 8 , and WFIRST 9 ) promise to provide 
precision cosmic shear measurements. 

Cosmic shear measurements are limited in precision by several systematic effects. It is important to understand and 
control these systematic effects in order to fully explore the potential of this probe (see for example Bacon et al. (2001); 
Bernstein & Jarvis (2002); Catelan, Kamionkowski & Blandford (2001); Croft & Metzler (2000); Erben et al. (2001); 
Heavens, Refregier & Heymans (2000); Heymans et al. (2004); Hirata & Seljak (2003a); Ishak et al. (2004); King & 
Schneider (2002); Refregier (2003); Takada & White (2004); Van Waerbeke & Mellier (2003); Brown et al. (2002) and 
references therein). One of the serious systematic effects of lensing is the correlated intrinsic alignment of galaxies which 
contaminates the lensing signal and acts as a nuisance factor (see for example Brown et al. (2002); Blazek, McQuinn & Seljak 
(2011); Catelan, Kamionkowski & Blandford (2001); Crittenden et al. (2001); Croft & Metzler (2000); Heymans & Heavens 
(2003); Hirata & Seljak (2003b); Jing (2002); Krause & Hirata (2011); Bridle & King (2007); Joachimi & Bridle (2010); 
Hirata et al. (2007); Hirata & Seljak (2004); Mandelbaum et al. (2006) ; Heymans et al. (2006); Faltenbacher et al. (2009); 
Okumura T., Jing (2009); Joachimi et al. (2010); Semboloni et al. (2008); King (2005); King & Schneider (2002, 2003); 
Kirk, Bridle & Schneider (2010) and references therein). For example, Bridle & King (2007); Joachimi & Bridle (2010) 
showed that if intrinsic alignment is ignored the determination of the dark energy equation of state is biased by as much as 
50%. Hirata et al. (2007) found that the matter power spectrum amplitude can be affected by intrinsic alignment by up to 
30%, showing the importance of developing methods to isolate the intrinsic alignment and remove it from the cosmic shear 
signal. 

There are two 2-point intrinsic alignment correlations. The first is a correlation between the intrinsic ellipticity of two 
galaxies, known as the II correlation. If the two galaxies are spatially close, they can be aligned by the tidal force field of 
the same nearby matter structure. The second intrinsic alignment correlation, known as the GI correlation, was identified by 
Hirata & Seljak (2004) and is due to a matter structure both causing the alignment of a nearby galaxy and contributing to 
the lensing signal of a background galaxy. This produces an anti-correlation between the cosmic shear and intrinsic ellipticity, 
since the tidal force and gravitational lensing tend to align the galaxy shapes in orthogonal directions. The GI correlation 
has been measured in various subsets of the SDSS spectroscopic and imaging samples by various groups. A detection of the 
large-scale GI correlation in the SDSS was reported by Mandelbaum et al. (2006) and then Hirata et al. (2007) found an 
even stronger GI correlation for Luminous Red Galaxies (LRGs). It was shown in these papers that this contamination can 
affect the lensing measurement and cosmology up to the 10% level and up to 30% in some cases for the matter fluctuation 
amplitude. This finding was confirmed by numerical simulations, where a level of contamination of 10% was found (Heymans 
et al. 2006). Further measurements of the GI correlation were made in the SDSS dataset by Faltenbacher et al. (2009); 
Okumura T., Jing (2009). Most recently, Joachimi et al. (2010) measured strong 2-point intrinsic alignment correlations in 
various SDSS and MegaZ-LRG samples. 

In a similar way, when we consider three galaxies and the related 3-point correlation, the cosmic shear signal (GGG 
bispectrum) also suffers from contamination by the 3-point intrinsic alignment correlations. The first is the III correlation 
between intrinsic ellipticities of three spatially close galaxies which are intrinsically aligned by a nearby matter structure. 
The second is the Gil correlation, where two spatially close galaxies are intrinsically aligned by a nearby matter structure 
which contributes to the lensing of a third galaxy in the background. Finally, there is the GGI correlation, where two galaxies 
are lensed by a structure which intrinsically aligns a third galaxy in the foreground. Unlike the 2-point correlations, the 
sign of the GGI and Gil correlations depend both on triangle shape and scale. Semboloni et al. (2008) showed that lensing 
bispectrum measurements are typically more strongly contaminated by intrinsic alignment compared to the lensing spectrum 
measurements, and that the contamination from the 3-point intrinsic alignment correlation can be as large as 15 — 20% 

1 http://www.cfht.hawaii.edu/Scicncc/CFHLS/ 

2 http:// www. darkencrgysurvey. org/ 

3 http://sci.esa.int/euclid/ 

4 http:/ / www.naoj.org/Projects/HSC/ 

5 http://www.stsci.edu/iist/ 

6 http://www.jwst.nasa.gov/ 

7 http://www.lsst.org/lsst/ 

8 http:/ /pan-starrs. ifa.hawaii.edu/ 

9 http://wfirst.gsfc.nasa.gov/ 
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compared to the GGG lensing signal. Finally, 3-point intrinsic alignment measurements are not only useful for constraining 
their contamination to 3-point lensing measurements, but are also useful for constraining models of intrinsic alignments and 
therefore constraining the contamination to all lensing measurements (including 2-point correlations) which will dominate the 
science cases of upcoming surveys. 

While the II and III intrinsic alignment correlations can be greatly reduced with photo-z's by using cross-spectra of 
galaxies in two different redshift bins (see for example Refregier (2003)) so that the galaxies are separated by large enough 
distances to assure that the tidal effect is weak, this does not work for the GI, GGI, and Gil correlations which happen 
between galaxies at different redshifts and large separations. The GI correlation and methods for its removal have been the 
topic of several recent scientific publications and we review these briefly. Initially, some first suggestions were discussed by 
Hirata & Seljak (2004). King (2005) extended the approach of template fitting by King & Schneider (2002, 2003) to 
include a treatment of the GI correlation. Bridle & King (2007); Joachimi & Bridle (2010) investigated the effects of the GI 
correlation on cosmological parameter constraints by assuming a model of the GI intrinsic alignment that is binned in redshift 
and angular frequency with some free parameters that are marginalised over. Kirk, Bridle & Schneider (2010) performed 
a cosmological constraint analysis where modelling of intrinsic alignment was included, showing a significant effect on the 
amplitude of matter fluctuations. Using a geometrical approach, Joachimi & Schneider (2008, 2009, 2010) proposed a nulling 
technique to remove the GI intrinsic alignment contribution by exploiting the redshift dependence of the correlations, but it is 
found that the technique throws out some of the valuable lensing signal. Most recently, the nulling technique has been applied 
at the 3-point level for the GGI correlation, but again with similar signal loss to that at the 2-point level (Shi, Joachimi & 
Schneider 2010). 

Finally, Zhang (2010a) proposed a technique to self-calibrate the GI intrinsic alignment signal by using the intrinsic 
galaxy ellipticity-galaxy density correlation, which requires that in addition to the galaxy ellipticity-ellipticity correlation 
(cosmic shear), one should also extract galaxy density-density and galaxy ellipticity-density correlations from the lensing 
survey. The GI correlation is then calculated and removed from the lensing signal. Most recently, Zhang (2010b) showed 
that redshift dependencies of intrinsic alignment can allow further improvements to the calculation of the intrinsic alignment 
contamination. The technique is commonly referred to as self-calibration because it uses correlations that can be extracted 
from the same gravitational lensing survey and used in order to calculate the GI contamination to the cosmic shear signal 
and remove it. Joachimi & Bridle (2010) applied an approach like the self-calibration, using correlations between lensing, 
intrinsic alignment, number density and magnification effects to constrain cosmological parameters. They found that the extra 
information from the additional correlations can make up for the additional free parameters in the intrinsic alignment so that 
the contamination can be removed without loss of constraining power. 

We organize the paper as follows. In Sec. 2, we briefly discuss the necessary survey parameters and lensing formalism. We 
then provide a summary of the 2-point GI self-calibration technique of Zhang with independent results. In Sec. 3, we develop 
the 3-point GGI self-calibration. We first establish an estimator to extract the galaxy density-density-intrinsic ellipticity 
correlation (ggl) from the observed galaxy ellipticity-density-density measurement for a photo-z galaxy sample. We then 
develop a relation between the GGI and ggl bispectra, which allows for the estimation and removal of the GGI intrinsic 
alignment correlation from the cosmic shear signal. Section 4 describes the residual sources of error to the GGI self-calibration 
technique, and we present the necessary relations to quantify these errors. These are compared to other sources of error in 
the bispectrum. Finally, we summarise the effectiveness and impact of the GGI self-calibration in Sec. 5. In the Appendix we 
expand upon the detailed calculation of the coefficients in the error calculation found in Sec. 4.1 and provide a list of typical 
expected values. 

2 BACKGROUND 

2.1 Survey information and weak lensing 

As mentioned in the previous section, the self-calibration technique proposed by Zhang (2010a) makes use of the information 
already found in a lensing survey (Bernstein 2009), including galaxy shape, angular position and photometric redshift, in 
order to calculate and remove the dominant intrinsic alignment contamination. In our performance calculations, we consider 
survey parameters to match a survey similar to those of the LSST lensing survey (LSST Science Collaborations and LSST 
Project 2009), but this is just an example. The GGI self-calibration technique is survey independent and can be applied to 
reduce the intrinsic alignment contamination in all lensing surveys. Galaxies are assumed to be sufficiently large and bright 
to be suitable for cosmic shear measurements, so we restrict any discussion of the self-calibration to these galaxies in order to 
avoid any sample bias. Galaxies are split into photo-z bins according to photo-z z p , where the i-th photo-z bin is described 
by a mean photo-z Zi and has a range Zi — Az;/2 ^ z p ^ Zi + Az;/2. In our notation, i < j implies that z; < Zj. The galaxy 
redshift distribution over the «-th redshift bin is nf (z p ) and rii(z) as a function of photo-z and true redshift, respectively, 
which are related by the photo-z probability distribution function p(z\z p ). 

In evaluating the performance of the self-calibration technique, we will consider as an example survey parameters to match 
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an LSST-like weak lensing survey (LSST Science Collaborations and LSST Project 2009), but of course the calculations are 
applicable to all current and planned weak lensing surveys (e.g. CFHTLS, DES, EUCLID, HSC, HST, JWST, LSST, Pan- 
STARRS, and WFIRST). We assume a survey coverage of half the sky (f s k y = 0.5) with a total galaxy surface density of 40 
and redshift density distribution of 



2 

Z 



with «o = 0.5. The ellipticity shape noise is described by 7 rms = 0.18 + 0.042z and photo-z error by a Gaussian probability 
distribution function (PDF) 



p(2|zP) = 7^T xp I V *W ' J' (2) 

with a z — 0.05(1 + z). We define photometric redshift bins of width Az = 0.2, centred at z; = 0.2(i + 1) (i — 1, • • • ,9). 
We do not include redshifts below z p = 0.3, not because of poor performance in the self-calibration, but rather due to the 
weaker lensing signal at lower redshifts. This artificially increases the fractional errors we evaluate in Sec. 4 with respect to 
the GGG lensing signal, as is evident in the increasing errors at low redshift in Fig. 3, and is not useful in evaluating the true 
performance of the self-calibration. 

The self-calibration technique relies upon two basic observables from a weak lensing survey. The first is the galaxy surface 
density of a photo-z bin, 5 E , which is a function of the 3D galaxy distribution S g . The second necessary observable is the 
galaxy shape, expressed in terms of ellipticity, which measures the cosmic shear 7. However, this cosmic shear signal is heavily 
contaminated by the intrinsic ellipticities of galaxies. There is a random component to this intrinsic ellipticity, which is simple 
to correct and which we ignore in the self-calibration calculations except as part of the shot noise in the error estimations of 
Sec. 4. A second component of the intrinsic ellipticity is due to the intrinsic alignment of galaxies caused by the gravitational 
tidal forces of large scale structure and was introduced in Sec. 1. 

We will label the measured shear as 7 s = 7 + 7 7 , where 7 7 denotes the correlated part of this intrinsic ellipticity due to 
the intrinsic alignment of galaxies. Since we are concerned only with the weak limit, we will work with the lensing convergence 
k instead. Thus from the measured 7 s , we obtain k 3 = k + k 1 . k is the projected matter over-density along the line of sight. 
For a fiat universe in the Born approximation, the convergence k of a source galaxy at redshift zg and direction 9 is 

«(*)=/ S(XL,0)W L (xL,XG)d X L. (3) 
Jo 

Wl( X l, Xg) is the lensing kernel and S(xl, 0) is the matter over-density in direction 8 and at comoving distance xl = X l{zl)- 
Xg = Xg{zg) is the comoving distance to the source. The comoving distance x is in units of c/Ho, where Ho is the current 
day Hubble constant. We assume a flat ACDM cosmology with Q m = 0.27 and S7a = 0.73. The lensing kernel is then 

W L (XL,XG) = ^ m (l + Z L ) X L(l-^) (4) 

2 XG 

when zl < zg and zero otherwise. 

In our calculations, we will work in Fourier (multipole I) space with the corresponding spectra to the correlations which 
can be built from these survey observables. The two-point correlation function is then related to the angular power spectrum 
and the 3-point correlation function to the angular bispectrum by 



<K(/i)S(/ a )> = (2n) 2 8 D (£ 1 +£ 2 )C(l 1 ) 
<ft(/i)S(/a)/s(/a)> = (27Y) 2 8 D (£ 1 +£ 2 +£ 3 )B(£ 1 ,£2,e 3 ), (5) 

where (• • • ) denotes the ensemble average and 8 D (£) is the Dirac delta function. For the bispectrum, 8 D (£i +£2 +£3) enforces 
the condition that the three vectors form a triangle in Fourier space. The 2D angular cross-correlation power spectrum is 
related to the 3D power spectrum and the 2D angular cross-correlation bispectrum to the 3D bispectrum through the Limber 
approximation 

B^M) - [ ^y' ^U;^, (6) 

where i,j,k denote the redshift bin and Wi(x) are weighting functions which depend on the quantity being correlated. For 
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example, when correlating weak lensing this is the weighted lensing kernel Wi(x) = Jo Wl{x\ X)fi(x')dx' where fi(x) is the 
comoving galaxy distribution in the i-th redshift bin. 



2.2 The 2-point GI self-calibration 

For the power spectrum, Zhang (2010a) proposed the GI self-calibration technique to calculate and remove the GI correlation, 
quantified by f 1 = C IG /C GG , from the angular cross-correlation power spectrum between galaxy ellipticity (k s ) in the i-th 
and j-th redshift bins. We confirm and summarise the important components of this technique here, as it is a necessary 
component in our own development of the 3-point GGI self-calibration. We also present independent performance estimates 
calculated using the linear alignment model for intrinsic alignment of Hirata & Seljak (2004) for the intrinsic alignment as 
described in Sec. 4. 

The following three observable correlations in a lensing survey between galaxy surface density and convergence for i < j 
are necessary to the GI self-calibration technique: 

cj?(i) « C GG (£) + C\f(t), 
C%\t) = C? G (l) + C°!(£), 

4 3) W = Cl?{l). (7) 

C*- 1 is the angular cross-correlation power spectrum between quantity a in the i-th redshift bin and fi in the j-ih redshift bin. 
a, /3 £ G,I,g, where G indicates gravitational lensing (k), I the correlated galaxy intrinsic alignment (k 1 ) and g the galaxy 
number density distribution (<5 E ) in the corresponding redshift bin. Thus we denote the GI power spectrum C' G in order to 
preserve the association of each quantity G or / to its redshift bin. By requiring i < j with sufficient photo-z accuracy, C GI 
and C'/ provide negligible contribution to the lensing signal and are neglected in C^K From these observables, Cl G (£) is 
expressed by using a deterministic galaxy bias b\(£) (see Appendix A of Zhang (2010a)) through the scaling relation 

Cl ^~ mk°' m (8) 

where W tj = J °° dxL J °° dxG[W L (XL,XG)fi{XL)fj{xa)] and IT, = J °° ff(x)dx- fi{x) = m(z)dz/dx is the true comoving 
distance distribution of galaxies in the i-th redshift bin and Wl(xl, Xg) is the lensing kernel. Using the effects of lensing 
geometry, C'f(l) is isolated from the second observable (£) (see Appendix B of Zhang (2010a)) using the estimator 

with Qiijl) = C? G \s(£)/Cf G (£). The subscript S denotes the correlation between only those pairs with zq < Zg. Q2(£) 
then measures the relative suppression of the gG signal due to the orientation dependence of the lensing geometry, where 
Ci G \s{£) < Cl G (£). For spectroscopic (true) redshifts, Cf G \s{£) = and the estimator C-f(£) is simply C^ ] \s{£) = C 9 J{£), 
since the gl correlation is independent of orientation. The measurement error (see Appendix C of Zhang (2010a)) in this 
estimator is 

+ ('+(r=w)^^°4 ,10) 

The fractional error on C\ G (and thus the residual statistical error in the measurement of C GG ) due to AC^ 3 is 

/ Wij \AClf 

A/lJ - \bi(£)mj c GG ■ (11 > 

This is also the threshold contamination f?j — flf resh at which the GI self-calibration will function at S/N=l. Similarly, the 
scaling relation in Eq. 8 is not exact, and its accuracy is quantified by 

sys _ bjjPjYLjj Cj G (£) 

ij W iS Clf{£) [ ] 

This introduces a residual systematic error in the measurement of C GG of dfij = e^J s f'j . 
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Figure 1. The inaccuracy of the scaling relationship in Eq. 8 is quantified in Eq. 12 by eJJ'^ . This inaccuracy is the source of the 
dominant systematic error in the measurement of C|? G due to the GI self-calibration technique. Left: e^ s is plotted for adjacent redshift 
bins, where the stronger dependence of the lensing kernel on redshift causes a significantly higher inaccuracy. Right: e^V s is plotted for 
redshift bins of varying distance from each other. As expected, the inaccuracy for these bin choices is generally less than for three adjacent 
choices. Using the linear alignment model for intrinsic alignment of Hirata & Seljak (2004), we find that our result is consistent with 
the toy model of Zhang (2010a) which gives the smallest systematic error. We thus expect a suppression of the GI intrinsic alignment 
contamination by at least a factor of 10 for adjacent bins and up to a factor of 50 for other bin pairs. These results are insensitive to the 
original GI contamination, such that for any /H 1 res " < fL < 1, the GI self-calibration will reduce the GI contamination down to survey 
limits or by a factor of 10 or greater, whichever is less. 



The third observable C\*>{t) can be used to calculate the galaxy bias. This gives a result for Cfj(£), which can then be 
removed from the first observable (£)■ Our notation above, which we will use throughout this paper, is slightly different 
from the original notation of Zhang (2010a) in order to be compatible with the GGI self-calibration. We denote analagous 
quantities in the 2- and 3-point self-calibration by the same variable, differentiated by the number of its indices. 

The GI self-calibration technique converts a systematic intrinsic alignment contamination /y into a residual statistical 
error A/ij < //, which is insensitive to the intrinsic alignment contamination. We find good agreement with Zhang's estimation 
of the performance of the GI self-calibration, based on independent calculations of the 2-point errors following the methods 
described in Sec. 4. Figure 1 shows the residual systematic (Sfij ) error in the measurement of C GG with the GI self-calibration 
technique. We use the linear alignment model for intrinsic alignment of Hirata & Seljak (2004) and find that our result is 
consistent with the toy model of Zhang (2010a) which gives the smallest systematic error. We thus expect for an LSST-like 
lensing survey a suppression of the GI intrinsic alignment contamination by at least a factor of 10 for adjacent bins and 
up to a factor of 50 for other bin pairs. These results are insensitive to the original GI contamination, such that for any 
^thresh ^ ji_ ^ ^ e q\ self-calibration will reduce the GI contamination down to survey limits or by a factor of 10 or 
greater, whichever is less. 
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3 3-POINT GGI SELF-CALIBRATION 

There are several sets of correlations between the observed galaxy surface density and convergence which can be constructed 
for galaxy triplets. Only three of these observed correlations are needed for the GGI self-calibration technique. The first is 
the angular cross-correlation bispectrum between galaxy ellipticity (n s ) in the i-th, j-th and fc-th redshift bin 



B^ k (e ly e 2 ,e 3 ) = B^ GG (£ U £ 2 ,£ 3 ) + Blf k G (l^£ 2 ,£ 3 ) + (2 perm.) + Bl ] I k G {£ 1 J 2 ,£ 3 ) + (2 perm.) + B$f(*i, £ 2 , t 3 ). (13) 

^tjk i s ^ ne angular cross-correlation bispectrum between quantity a in the j-th redshift bin, /3 in the j-th redshift bin and 7 in 
the k-th redshift bin. a,/3,j £ G,I,g, where G indicates gravitational lcnsing (k), / the correlated galaxy intrinsic alignment 
(k 1 ) and g the galaxy number density distribution ((5 E ) in the corresponding redshift bin. Unless catastrophic photo-z errors 
overwhelm the data, we can safely neglect the correlations Gil and III which require spatially close galaxies by selecting 
galaxy triplets where i < j < k. Under this requirement, we also have B/, GG 3> B GGI , B GI k G due to the lensing geometry. We 
then have for i < j < k, 



B^ k (£i,£ 2 ,£ s ) w B ggg (£i,£ 2 ,£ 3 ) + B I i G k G (£ 1 ,£ 2 ,£ 3 ). (14) 

Thus the dominant intrinsic alignment contamination is from the GGI bispectrum B\ G k G (l\,l 2 ,£ 3 ) (i < j < k), which the 
GGI self-calibration technique seeks to calculate and remove. Here we denote the GGI bispectrum B IG k G in order to preserve 
the association of each quantity G or / to its redshift bin, as for the 2-point GI power spectrum. 

The second correlation is measured in the angular cross-correlation bispectrum between convergence (k s ) in the j-th 
redshift bin and galaxy density (5 E ) in the j-th and k-th redshift bins. Of interest to the self-calibration is the case where 
i = j = k, and we have 



B™(t u t2,t 3 ) = B^itut^ + B^ (t u t a ,£ 3 ). (15) 

This correlation contributes further information about the intrinsic alignment of galaxies 

The final correlation of interest is measured in the angular cross-correlation bispectrum between galaxy density (<5 E ) in 
the i-th, j-th and fc-th redshift bins when i = j = k, giving 



Bi:>(£ 1 J 2 ,£ 3 ) = B%°(£ 1 ,£ 2 ,£ 3 ). (16) 

We also require for the GGI self-calibration those observables in Eq. 7 for the GI self-calibration. It is important to 
note that we have thus far neglected the contribution of magnification bias to these measurements. This will be further 
discussed and justified for the 3-point measurements in Sec. 4.3 and was discussed and shown to be negligible for the GI 
self-calibration by Zhang (2010a). There is also a non-Gaussian contribution to the observed bispectra. We briefly discuss 
the impact of this non-Gaussianity on the self-calibration technique in Sec. 4.4, but otherwise leave discussion and calculation 
of this non-Gaussian contribution to the bispectrum to other works. 

Our GGI self-calibration technique will calculate and remove the GGI contamination in Eq. 14 by using the measurements 
from Eqs. 15 & 16, which are both available in the same lensing survey. We express the fractional contamination to the lensing 
signal by the correlated intrinsic alignment as 



f i (f f „ n _ B I i G k G (£ 1 ,£ 2 ,£ 3 ) . . 

fijk(£i,£2,£ 3 ) = B^ GG (£x,£ 2 ,£ 3 y (17) 

For the self-calibration to work, the contamination flj k (£i,£ 2 ,£ 3 ) must be sufficiently large as to contribute a detectable 
B]*?{£ u l 2l £ 3 ) at the corresponding £ bins in B$ (£ u £ 2 , £ 3 ). We denote this threshold f$ k reah . When f? jk > f^ k resh , the GGI 
self-calibration can be applied to reduce the GGI contamination. The residual error after the GGI self-calibration will be 
expressed as a residual fractional error on the lensing measurement. In our notation, we differentiate Afij k as statistical error 
and Sfijk as systematic error. The performance of the GGI self-calibration will then be quantified by the parameters fij" k eah , 
Afijk and Sfijk, which are discussed and calculated in Sec. 4. 



3.1 Relationship between B\ G k G and B\?? 

The first step in the GGI self-calibration is to determine the relationship between B{f k G and B?? 9 . Under the Limber approx- 
imation, the 2D GGI angular cross-correlation bispectrum between the j-th, j-th and k-th redshift bins is related to the 3D 
matter-matter-galaxy intrinsic alignment bispectrum by 
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Jo X \ X X ^ 



where 



Wi( X L)= W L ( X L,XG)MXG)dXG. (19) 

The integral runs from zero to oo in order to take into account the photo-z error. We again denote the GGI bispectrum B\f^ 
in order to preserve the association of each quantity G or I to its redshift bin and will continue this convention throughout the 
paper. Similarly, the 2D ggl angular auto-correlation bispectrum is related to the 3D galaxy-galaxy-galaxy intrinsic alignment 
bispectrum by 

B^(e 1 ,t a ,ts) = f ° O ^b 33iI (k! = pfc = |,*3 = |;x) d X . (20) 

We will adopt a deterministic galaxy bias b 3t k (Fry & Gaztanaga 1993) such that the smoothed galaxy density is a function 
of matter density expressed as 

5 g (x; X ) = & fl ,i(x)M*;x) + b -^^S 2 m ( X ;x) + 0(8 3 ). (21) 

The first term b g ,i is the linear galaxy bias (as used by Zhang (2010a) for the 2-point correlations). The second term represents 
the first order non-linear contribution. b 3i2 is typically found to be negative and ^ 6 g ,i (Cooray & Sheth 2002). Unlike in 
the 2-point case, it is insufficient to model the bias as simply scale dependent (Jeong & Komatsu 2009). Following the 
galaxy-galaxy-galaxy halo bispectrum derivation of Jeong & Komatsu (2009), we use this expression of the galaxy density 
to relate B' s f^ to -B^y • We neglect the portion of the bispectrum due to primordial non-Gaussianity and the trispectrum 
term, which contains further information about the non-Gaussianity. This is justified and discussed further in Sec. 4.4. This 
results in the relationship 

B ggl i(k 1 ,k 2 ,k 3 ;x) = &g,i(x)-B<s<s 7 ' (fci, k 2 , k 3 ; x) 

+b g ,i(x)b g ,2(x) [P Sl i(ki;x)P6s(k2;x) + P66(k2;x)P Sl i(k 3 ;x) + P Sl i(ki;x)P s ^(k3;x)] • (22) 

If the galaxy bias changes slowly over the i-th redshift bin with median comoving distance Xiy we can write to a good 
approximation b' k — b 9t k{Xi)- Substituting Eq. 22 into Eq. 20, we have 

^f(£r,£ 2 ,4) = f°° %^((b\) 2 B SSlI (k u k 2 ,k 3 ; X ) 
Jo X v 

+b\bi [P Sjl (fci; x)P M (fe>; x) + Pss(k 2 ; x)^V (*s; x) + P«y (fci; x)P^i (fcs; x)] )d X - (23) 

We can further approximate B(ki, k 2 , k 3 ; x) ~ B(ki, k 2 , k 3 ; Xi) an d P{k;x) ~ P{k\Xi) m t- ne limit where the comoving 
distance distribution of galaxies in the i-th redshift bin is narrow. This leads to the following approximations of Eqs. 18 & 23, 

4f(<i,4<s)«V(' i i>fe.*»i*)% ( 24 ) 

Xi 

and 

Xi 

+b\bl [P Sj i (fci ; X i)Pss{k2; Xi) + ft«(fe»; X*)^V (*»; X<) + (*i; xO^V (*s; X*)] ) , (25) 
where W^ fc = / °° /i(x)W^'(x Wx)dx and IT;;; = / °° ff( X )d X . From Eqs. 24 & 25, we have 



\ u l) LL m u lXi 

x [P Sl '(ki;Xi)Pss{k2;Xi) + Pss{k2\Xi)P Sl '(k 3 ;xi) + P6 7 '(ki;Xi)P Sl '(k3;Xi)] ■ (26) 

In order to express the 3D power spectra in Eq. 26 as 2D spectra, we will use the approximation made by Zhang 
C'f(£) « P Sl i (fc; Xi) and the similar approximation C% G (£) w P^fc; X<)^i where cj^ = J °° Wj 2 (x)^X- Equation 26 is 
then 
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i ijk (luiiJs) ~ 7^-^ (4,^4)- 



Big- (t u t2, h) « TiT^-^ (^i.^.^3) 

— ^, T /.]_■_■ 

(27) 



<#(*i)Cg G (* 2 ) + Cg G (* 2 )C#(*s) + ^-C^C^h) 

This relationship, while developed in the same way as for the GI self-calibration, is necessarily more complicated due to 
the inclusion of the non-linear galaxy bias and the presence of the GG correlation. Thus in order to apply this relationship, it 
is necessary to not only develop an estimator for B^ 3 (ti, £2, £3), which we describe in Sec. 3.2, but also to use the estimator 
C*?(£) in Eq. 9 developed for the GI self-calibration and the resulting C GG (£), as measured by the GI self-calibration (Zhang 
2010a). The GGI self-calibration technique is thus dependent upon the resulting measurements of the GI self-calibration 
technique. 

3.2 B^ 99 Measurement 

Information about the galaxy density-density-intrinsic ellipticity bispectrum, B\ 99 , is contained within the observable B^f- — 
B\f? + B u 99 . To measure it directly, we must first remove the contamination of B^ 99 . For a spectroscopic galaxy sample, 
lensing geometry requires eliminating those triplets of galaxies where the redshift of the galaxy used to measure the ellipticity 
is lower than those used to measure galaxy number density. In this way, those triplets remaining have no contamination from 
B %! 9 and measure only B- 99 . 

In the case of a photo-z galaxy sample, this is not possible due to typically large photo-z error. Even for a photo-z bin 
with Az — > 0, the photo-z error causes a true redshift distribution of width 2o\p = 0.1(1 + 2). In practice, photo-z bins are 
typically > 0.2. With such large errors, it is possible for galaxy triplets in the i-th redshift bin to provide a measureable lensing 
contribution to -B G f 9 even when requiring that the redshift of the galaxy used to measure the ellipticity is lower than those 
used to measure galaxy number density, except for the special cases where we limit to sufficiently low values the redshift or 
both the photo-z error and bin size. A more careful approach is thus required when separating B{?? from £> G f 9 for a general 
photo-z galaxy sample. 

We apply the approach used by Zhang for the power spectrum to the bispectrum B^- , wherein we consider the 
orientation dependence of the two components. We will first define a redshift for each galaxy in the triplet: z G /i for the 
galaxy used in the lensing/intrinsic alignment measurement and z B , z g i for the two galaxies used in the number density 
measurement. The ggl correlation is independent of the relative position of the three galaxies. For example, the correlations 
with zi < z g < z g i , z g < zi < Z g i or z g < Z g i < zi are st at ist ically identical when the sides of the triangle are fixed. However, 
the Ggg correlation does depend on the relative position of the three galaxies. Due to the lensing geometry dependence, the 
correlation with zg < z g , z g i is statistically smaller than other orientations. 

This dependence provides two observables from B^ . The first is B$ , where all triplets are weighted equally. The second 
is B^-\s, which counts only those triplets with zg < z g , z g i. This weighting is denoted by the subscript 'S'. From our previous 
discussion, we then have B\ 99 = B'? 9 \s and -B G f 9 > Bf i 99 \s- We now define the ratio 



Q*(tut*,t*) = ( 28 ) 



B%°\s{£i,e2,£ 3 ) 
B°? 9 (£u£2,t 3 ) 

where we have explicitly included the ^-dependence which had been neglected previously in this section. This ratio describes 
the suppression of the signal due to the weighting of triplets described previously. By definition < Q3 < 1, with Q3 = if 
the photo-z is perfectly accurate and Q3 — 1 if the photo-z has no correlation to the true redshift. Q3 is calculated using the 
galaxy redshift distribution, which is discussed in Sec. 3.3. 

We now define an estimator for B^ 99 (that we denote -B^ 99 ) in terms of Q3(£i,£2,£ 3 ) and the two observables 

r(2)/ 



^,£ 3 ) + B^ 9 (tut2,£ 3 ), 

B™\ 8 {l u t2,h) = B{ 9 ?{tut2,£3) + Bf 99 \s(tu£2,£3). (29) 
This estimator is 



r,I 99(f f gj I S (ii ,£2 ,£3 ) ~ Q3 (h ,£2 , £ 3 )B% (h ,£ 2 ,l3) . . 

Biii (£l ' £2j:i) = R#M) • (30) 

As expected, when Q3 — this gives B\ 99 — B^j \ s as for a spectroscopic galaxy sample with no photo-z error. However, 
Q3 must not approach unity, where B\ 99 is singular. For the LSST-like survey described in Sec. 2.1, we calculate Q3 for 
various redshift bins following the procedure described in Sec. 3.3. This result is given in Fig. 2 for equilateral triangles, where 
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Figure 2. The behavior of Qz{h,(-2,h) = B°f 9 \s(£i ,£2, h)/B^f 9 (li , l 2 , £3) for equilateral triangles (£ = l x = £ 2 = £3) over three 
redshift bins spanning the survey range. Similar to the 2-point case, the suppression is dependent on the redshift bin chosen, increasing 
with redshift due to increased photo-z error at higher redshift, but is largely scale independent due to being the ratio of two bispectra. 
For this reason there is also little dependence on triangle shape. Generally, Q 0.4, and the significant deviation from unity ensures 
that the estimator B^ 99 is valid for lensing surveys of interest. 

we find Q3 w 0.4 and in general that Q3 should deviate significantly from unity. The estimator B\ 99 is thus expected to be 
applicable in any typical lensing survey. 

3.3 Evaluating Q 3 (ti,£2, £3) 

In order to evaluate the ratio Q3 in Eq. 28, we will begin from the real space angular correlation function w G " (61, 62, 63; Zq, Zg , 1 
between shear at photo-z zq and galaxy density at photo-z Zg and zp . The average correlation over the distribution of galaxies 
in the i-th redshift bin is 



'•"• < ""' ^O 1 ,e 2 ,03;zc,Zg,Zg,) n[ (zc)n[ {zg)n[ (z^,) 



dzQ I dzf I dz P 



ir';:r''uh.<i..ih) - / </:/; / d< j dz g :w L 

f /*r*n r-oo roc 

.„ , ....„■ , dzc / dzg I dZg,w G " (81,82,63; ZG,z g ,z gl ) 

Ji JO JO JO 

xp(zG\z G ')p(zg\zg)p(zg,\zg l )nf'(z G ')n P (Zg)nf'(zg,), (31) 

where we have used the shorthand f. — fg'^&^M ^° represent integration over the i-th redshift bin. In terms of the ensemble 
average (•••}, which is in practice an average over 8' , the angular real space correlation function is 



poo 

° gg (81, 82,83; ZG,Zg,z g ,)= / dz L (S m (8';zG)S g (8' + 8 2 ;zg)Sg(8' + 8 3 ;Zg,))WL(zL,ZG)S D (8 1 + 8 2 + 83), 



where S D (8 1 +82 + 83) ensures that the three vectors form a triangle. We can now write Eq. 31 as 



(32) 



f f f poo roo roo r oo 

w ui 9 (81,82,83) = dzc dzg I dz P t I dz G I dz g I dz gl / dz L (8 m (d' ; z G )S g (d' + 8 2 ; z g )S g (d' + 83; z gl )) 

Ji Ji Ji JO JO JO JO 

xS D (8i + 8 2 + 83)W L (zL,ZG)p(zG\zG)p(zg\zg)p(zg l \zP)n[(zG)ni(Zg)nf'(Zg l ) 



dz L / dz g / dzg,(5 m (8';z G )Sg(8' + 8 2 ;Zg)Sg(8' + 83;z g ,))S D (8i+ 8 2 + 83) 



xWi(z L )m(zg)ni(z gl ). 



(33) 
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The second correlation function needed is identical to Eq. 33, but takes the average over all triplets such that z G < z p , z p ,, 



poo p oo poo 

™?u g \s(0i,02,e 3 ) = / dz L dzg dz g ,(8 m (e'; ZG )Sg(e' + e2;z g )s g (8' + e 3 ;zg,))s D (e 1 + e 2 + 8 3 ) 

Jo Jo Jo 

xWi(z L )ni(zg)ni(Zg,)rj(z L , z g , z g >). (34) 

We have used here 

, . _ 3 / t dz G Sj dz g Sj dz g> So°° dz G W L (zL, z G )p(z G \zG)p{z g \Zg")p{zg.\Zg,)nf (zg)raf \z P )nf (z p ,)S(zg, z P , z p ,) 

V(ZL ' Z3 ' Z3 ' ] ~ Si dz a Si dzf Si K fo°° ^^(zi.zoJp^lzgJp^l^fJpCvl^Jnf (^S)nf (^f)nf (<,) 

(35) 

where S(z G , z g , z gl ) = 1 if z G < z g , z p , and is zero otherwise. Since S(z G , z g , z p ,) allows only 1/3 of the integral to survive, 
t)(zl, z g , z g r) is normalised by a factor 3 in order to remove the suppression due to the selection function and measure only 
that due to the lensing geometry. This is demonstrated by the relation 

J. dz p G f. dz p g f. dz p g , /„" dz G p{Z G \z^)p{Z g \z p )p{Z g ,\z p )n p (z^)n p {z P )n p (z p )S(z^, zj, z p g ,) 1 

f.dz* Sidz P Sidz p ,S~ dzGp(Z G \z p )p(Z g \z p )p(Z g ,\z p )n p (z p )n p (z P )n p (z p ^ 3' ( > 

We now take the Fourier transform of Eqs. 33 & 34 to find the bispectra Bfj 9 and Bf i f 9 \s, respectively. Again following 
the Limber approximation, with dominant correlation at zl — z g = z g i, we have 



r B^ {kl MM-,x) m{X)f ; {X) dx (37) 
Jo X 



and 



B% a \s (li , £2, is) = r B G " (k, ,k 2 ,k 3 ; X ) Wi{X }J\ {X \ { X , x(z g ) = X,x(z 3 >) = x)dX- (38) 
Jo X \z) 

The ratio Q 3 is now expressed directly through Eqs. 37 & 38. We can approximate Q 3 w fji, where fji is the mean value 

of 77 across the «-th redshift bin, since the integrals differ only by a factor 77. r\ has the same dependence as Q 3 on the relative 

contribution to the Ggg correlation from triplets with z G < z g , z p , compared to triplets with other relative orientations. In the 

limit where photo-z error dominates, o P 3> Az, and there is no suppression of the contribution to the Ggg correlation by the 

selection function, so 77, Q 3 — ¥ 1. In this limit, the estimator B\f 9 becomes singular and B^ 9 can no longer be differentiated 

from B^ 99 . In the opposite limit, where a p <g; Az, the selection function suppresses all contribution to the Ggg correlation 

and 77, Q 3 — > 0, where our estimator mirrors the extraction method for B\f? in spectroscopic galaxy samples. 



4 PERFORMANCE OF THE GGI SELF-CALIBRATION 

In order to evaluate the statistical and systematic errors in the GGI self-calibration, we calculate directly the power spectra 
and bispectra through the Limber approximation according to the anticipated survey parameters discussed in Sec. 2.1. For 
the bispectra, we employ the fitting formula of Scoccimarro & Couchman (2001) for the 3D matter density bispectrum due 
to non-linear clustering. We modify this as described in Sec. 3.1 for the 3D galaxy bispectrum, using values for the galaxy 
bias of 61 = 1.0 and b\ = —0.1 (Simpson et al. 2011). We include the intrinsic alignment correlations of B IGG and B 1 " in 
a straightforward manner following the linear alignment model of Hirata & Scljak (2004), where Pg^i = — D (z)(i+z) 
Bridle & King (2007), we extend this to the non-linear matter power spectrum for use in the fitting formula, where C\ is 
estimated by comparison to Fig. 2 of Hirata & Seljak (2004). 



4.1 The estimation of B] 99 

In order to quantify the accuracy of the estimator B\ 99 {1\ , £2, £ 3 ), we consider the contribution of measurement errors such as 
shot and shape noise in B^(£i, £2, £3) which propagate into our measurement of B\??(£\, £2, £ 3 ) through the estimator. We 
calculate the rms error for a given redshift bin, working in a pixel space with Np sufficiently fine and uniform pixels of photo-z 
with bin width Az and angular position with bin width A£. Each pixel is associated with a photo-z z p , angular position 9 a , 
measured overdensity 5 a + an d measured 'shear' K a + k 1 + , where 'N' represents the measurement noise. From Eq. 29, 
we construct the pixel space angular bispectra 
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S (2) (4,£ 2 ,£ 3 ) = Np 3 Y,[S a + 5%][S + S^}^ + 4 + k?] exp[t(*i ■ 6> a + ^ 2 ■ 9 P + £ 3 ■ 7 )], 

a/37 

B (2) |s(<M2,4) = JVp'^Ja+^l^+^l^+K^+K^lexp^l •ea+/a-e/9+/8-eT)]Safl 7 . (39) 

ck/3 7 

S a ^ 7 = 1 when , Zp > 2 7 and is zero otherwise. Thus in the limit Np S> 1, $^ a(97 £ a/ 3 7 = Np/3 and the average S a( s 7 = 1/3. 
From our definition of the estimator in Eq. 30, we can construct the difference 



- B lu = (1 \ s) N p 3 E ex P[^ -Oa+ta-Op + ia- 0^][{6 a + 5^)(5 P + 8f)fa + k\ + ^)(3S a ^ - Q 3 ) 

a/3 7 

— (1 — Qz^aSpulf] 

= jz — -p^rrNp 3 \^ exp[i(£i ■ a + £2 ■ Op + £3 ■ 0~t)][{5 a 5p + S^Sp + 5™ <5^)(k 7 + k 7 + k%) + 5 a 5p(Ky + k%)] 
x(3S a/37 -Q 3 ). (40) 

Here we have used S a py = 1/3 and that the ggl correlation doesn't depend on the relative position of the galaxy triplets. The 
rms error is 



(AB't? 9 Y = 1 Np 6 E ex P[^i • 0a + *a • 0/a + *s ■ 7 )] eap[i(/i ■ X + £2 ■ 6>„ + £3 ■ 0„)}(3S aM - Q 3 ) 

apy A^ 

X(3Sa^ - Q3)([(5 a 5^ + &a 6p + 5^ <5^)(K 7 + K 7 + K^) + 5 a 5^(K 7 + Ky)] 

x [(5a<5^ + 8\5 V + 8^8^)(k v + k v + k") + 5\8y,(Kv + k„)]), (41) 

where (• • • ) is the ensemble average. The ensemble average is over 121 terms of the form (ABCDEF), A,B,C,D,E,F G 
To simplify this we apply Wick's theorem for the 6-point correlation, 



(ABCDEF) = (AB) (CD) (EF) + (AB) {CE) (DF) + (14 perm.) . (42) 

This results in 1815 products of three 2-point correlations, most of which are zero. Any correlation between signal and 
noise or dissimilar noise terms vanish. Due to the angular dependence of the correlations ((A a Bb) = WAs(#a — Ob)), only those 
correlations with (A a Bt) where a € a, /?, 7 and b £ \,fi,v are non-vanishing. This leaves 42 surviving products: 



(aB^Y = l — 2 N P 6 V V exp[t(*i ■ a + £2 ■ Op + £3 ■ fty)] eap[t(/i ■ X + £2 ■ 0„ + £3 ■ 0„)] (35 a(37 - Q 3 ) 

x (35a m „ - Q3) [(<Wa) [((fytf/j) + (<« 7 «i/) + ( K 7 / ' t ^)) + (^^)(« 7 «i) + («7*/»)(<5j9«i/)] 

+(<5 a <5 M ) [((5/s5a> + (SpSx)^ + (KyK^)^j + (8p8x) (n^nl) + (n~ f 6x)(6pK v )^ 

+($a Sx) [((5/3<5/»> + {Sp8^)^ ([k~,k v ) + (^k^) + (k 7 k£)) + ((k 7 + J 7 ) <5 M )(^ + 

+ (<5a <5^) [((<5/3<5a) + {SpSx)^ ({KyKv) + (HyK^) + (K 7 K^)^ + (^K 7 + 5a)(<5/3 (Ki/ + A/))] 

+(5a^) [(<5^<5 m )(k 7 <5a) + (<5^<5a)(k 7< 5 m )] + (<5 a + n^j) [(Sp 5^){^, + n 1 ^ 8 X ) + (<5^<5 A vr )(^K 7 + n 1 ^ <y|ij) 

Noises only correlate at zero lag ({5^ 5 X ) oc 8 a \, (k™ k^) oc 5-y V ), and the correlations (55), {8k 1 ), (kk) and (k 1 k 1 ) depend 
only on separation, not on relative orientation of the galaxy pairs along the line-of-sight. However, (kS) is dependent on the 
relative orientation along the line-of-sight and must be treated with care when evaluating Eq. 43. In order to quantify this 
orientation dependence, we apply Q2 such that 

(W ^K(T^) + lf) (W - (44) 

We can now evaluate Eq. 43 analytically, taking the Fourier transform to find 
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Figure 3. Left: The residual statistical uncertainty A/W in the measurement and threshold of intrinsic alignment contamination 

fl^ esh at which the GGI self-calibration technique can calculate and remove the intrinsic alignment contamination at S/N=l are plotted 
for a variety of redshift bin combinations. Right: The minimum measurement error for B^GG ; s plotted for comparison. Both errors 

are plotted for equilateral triangles (£ = £1 = £2 = £3) and have a similar i dependence, with the effects of shot noise taking over at large 
£. Generally, A/^ < an d is thus negligible. We expect this result to hold for non-equilateral triangles as well, but the use of the 

GGI self-calibration is limited by our understanding of non-Gaussian effects for very elongated triangle shapes, as discussed in Sec. 4.4, 
and we leave discussion of its applicability for these very elongated triangle shapes to a future work. 



(A^f) = 2{(c GG C°? + bC GG ' N C°? +2fCfcf)c?f 

+2(a [(c% G + cl!) c?? + c?*c?f\ + dc GG > N c°? + g cfcf + ihcfcU 

+ [c (c GG + CU) + eC GG ' N ] r;r V )c,"/" V }. (45) 

The details of this calculation and the coefficients a — h are found in the Appendix. 

The final rms error AB^ff evaluated for a given triangle with bin width A£ is then given by 



/ar 7 99^ 2 — — / fr< GG r<s3 I ur< Ga ' N nao i o fn9 G n9 G \ n gg 

\r Uiii ) ~ £ 1 £ 2 e 3 Ae 1 Ae 2 A£ 3 f sky 1 °- + bC - Gii+i/LiiL »j Gii 

+2(a [(C GG + cU)c% + CgC$] + dC GG ' N Cff + gCfcf + 2hCf Of* 

+ [c (C GG + CU) + eC GG - N ] C!t N )c!t N ). (46) 

C 99,N = l/n . C GG,N = 7 3 M / n . j where 

rii is the average number density of galaxies in the i-th redshift bin. Unlike the 
GI self-calibration, AB^ff is dependent on the intrinsic alignment contamination through C\f . However, it is still insensitive 
to the intrinsic alignment contamination in the limit where C GG is dominant. 
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The errors AB\f? and AC\? propagate into the measurement of B\ GG through Eq. 27. Performing a standard error 
propagation gives a residual statistical error AB{ GG . For the equilateral case, this simplifies to 




where we have neglected terms of order A 2 . To find the fractional error Af^ k this induces in the lensing bispectrum, we 
simply scale ABl G k G by the factor fl jk such that Af^ k = f\ ik AB\ GG . Like the 2-point case, this error is equal to fl^ k esh , 
the minimum intrinsic alignment flj k which can be detected through the self-calibration with S/N=f or AB^ 99 = B\?? '. Thus 
fijk esh = A/ 4 " fc represents for the self-calibration both the residual statistical error in the measurement of C GGG and the 
lower limit at which the intrinsic alignment can be calculated and removed. The GGI self-calibration technique can then turn 
a systematic contamination f' jk of the lensing signal into a statistical error Af^ k < f' jk which is insensitive to the original 
intrinsic alignment contamination. 

We compare the error Af.^ to the minimum rms error due to cosmic variance and shot noise in the B GGG measurement, 
which ignores other sources of error like the intrinsic alignment. The rms error of B GGG (i / j / k) is 



( AD GGG\ 2 ( n GG . n GG,N\ ( n GG . n GG,N\ ( n GG , n GG,N\ . I ' n GG , n GG,N\ n GG n GG 

y±&ijk J - [yu + <^u J [yjj + j [ykk + ^ kk j + [yu + <~>u j ^jk <^jk 

i / j^iGG | j^tGG ,N \ *~<GG f-iGG , / ^<GG , ^iGG,N \ *~<GG ^iGG , c\/~iGG /~<GG *~iGG / aq\ 

+ I ^33 +L/ jj j^ik l^ik +[<^kk + ^ kk j L 'i3 + ° jk <~'ik ■ 

This gives an absolute lower limit on the fractional measurement error of 



' eyfc ) WsA/i^Aft/rt, B GGG 1 ■ 



ijk 

Where Af^ k < tijk, the residual measurement error introduced after the GGI self-calibration is negligible, with very 
little loss of cosmological information. We find this to be true for an LSST-like survey, as shown in Fig. 3. More generally, 
since Af^ k and etj k scale similarly with respect to survey parameters, this should hold for other lensing surveys as well. 



4.2 The accuracy of the Bffi-B-? 9 relation 

In addition to the measurement error introduced through the estimator B{? 9 , there is a systematic error which is introduced 
by Eq. 27, which relates the intrinsic alignment contamination B\f k G in the lensing bispectrum to other survey observables. 
The accuracy of Eq. 27 is quantified by 



W ijk Bl^{£ u £2,e 3 ) b\ W ijk 



(&i) 2 n iii s/° G (£i,£ 2 ,4) (b\) 2 w«n« B I i GG (e 1 ,e 2 ,e 3 ) 

C\?{^)C GG {1*) + C GG (£ 2 )Ci 3 (£ 3 ) + ^C#(*)<#(*3) 



(50) 



This induces a residual systematic error in the lensing measurement of 



Sfijk = etjkfijk- ( 51 ) 
e tjk i s evaluated numerically and shown in Fig. 4 for equilateral triangles. As in the 2-point case, Eq. 27 is most accurate 
for those galaxy triplets which do not share neighbouring redshift bins. In the cases of neighboring bins, the lensing kernel 
varies more quickly due to the proximity of the galaxies in redshift. This causes Eq. 27 to be less accurate, increasing the 
systematic error. For galaxy triplets with bins which are not adjacent, |e^f | < 0.1. For these bin choices, the intrinsic alignment 
contamination can be suppressed by a factor of 10 or greater. In most cases where two or three bins are adjacent, \e*j£\ < 0.2, 
which allows for a suppression in the contamination by a factor of 5-10. In only a few of the cases where all three bins are 
adjacent is |e|j^| > 0.2, and even in these cases we expect a suppression in the contamination by a factor of 3 or more. 
These results are insensitive to the original intrinsic alignment contamination, such that for any flf k esh < f'j k < 1, the GGI 
self-calibration will reduce the GGI contamination down to survey limits or by a factor of 5-10 or greater, whichever is less, 
for all but a few redshift bin triplets. 
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Figure 4. The inaccuracy of the relationship between Bj^jp and the observable B^ 9 is quantified in Eq. 50 by e"V^. This inaccuracy is 
the source of the dominant systematic error in the measurement of B^ GG c j ue £ Q ^ ne qqj self-calibration technique. Left: e^Jjf is plotted 
for three adjacent redshift bins, where the stronger dependence of the lensing kernel on redshift causes a significantly higher inaccuracy. 
Right: e^Jjf is plotted for redshift bins of varying distance from each other. As expected, the inaccuracy for these bin choices is generally 
less than for three adjacent bins. e"Jfj* is plotted for equilateral triangles {I = i\ = £2 = ^3) in all cases. Equation 27 is usually accurate 
to within 20%, except for some adjacent bin choices, where it reaches a maximum of approx. 35%. Despite the inaccuracy of Eq. 27 
being greater than for Eq. 8 in the GI self-calibration, the GGI self-calibration is still expected to reduce the GGI intrinsic alignment 
contamination by a factor of 5-10 or more for all but a few adjacent redshift bin triplets. These results are insensitive to the original 
intrinsic alignment contamination, such that for any ff-£ e < ffj k < 1, the GGI self-calibration will reduce the GGI contamination 
down to survey limits or by a factor of 5-10 or greater, whichever is less, for all but a few adjacent redshift bin triplets. 



4.3 The magnification bias 

In addition to distorting the shapes of galaxies, gravitational lensing introduces a magnification bias to the observed galaxy 

overdensity Sg — S g + 2(a — l)/c, where a is determined by the logarithmic slope of the unlensed galaxy luminosity function. 

(21 

The magnification bias affects all three observable bispectra, but we expect the dominant contribution to occur in B^J . 
Including the average magnification bias in the i-th redshift bin, m* = (2(a — 1)}, Eq. 15 is modified to be 



4 s ? = + Blu + 2m, (B^ + B%') + m? (sgf ° + B^ G ) . (52) 

We seek to measure B^ff for the GGI self-calibration, so we will examine the effect magnification bias has on this 
measurement. Applying the estimator in Eq. 28, the bispectra Bjjj? and B GGG are unaffected, while the others are suppressed 
by a factor similar to (I — Q3). Thus the estimator acts to measure the dominant combination B\f? -t-m^ B G f G , where m%B GGG 
contaminates the B\f® measurement. Because the GGI self-calibration depends on the results of the 2-point self-calibration, 
we will also require the contribution Cff + m,iC GG from Cjf as discussed by Zhang (2010a). 

We cannot remove this contamination with any certainty due to measurement errors on m,i, Cfi and Bjfj . The direct 
estimation of the errors involved is lengthy, so we will instead determine the accuracy to which these measurements must 
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be made in order for the contribution due to magnification bias to be negligible with respect to other errors in the GGI 
self-calibration. 

We will assume m» has some measurement error Am;, B GGG a measurement error AB GGG and C GG a measurement 
error AC GG . From Eqs. 27 & 8, the induced measurement error in B\ G h G is 



AB^ = J|- (2 mj Am,B«f G + G A^ G ) " |^ + l) (Am,Cf + mi Cf ACf ) 2 . (53) 

Since we are only interested in the upper limit of this effect, we note that B GG > B GG for i < j < k and use the reduced 
bispectrum for equilateral triangles B GGG = 3QkC GG C GG to write a simplified expression for the induced fractional error in 
the B GGG measurement as 



< 



W; 



ijk 



(b\) 2 n v 



< 0(10" 



2 |m,Ami| + rrii 
Ami 



AB 



GGG 



B GGG 



b%W„ 



3Q k (b\rili V w, 



2 w i3 



+ 1 



[\Ami\ + 



ACh 



GG 



c GG 



0.1 



+ rrii 



a tdGGG / nGGG 
^^iii I J-'iii 



10% 



+ 



Ami 
0.1 



AC£; G /C^ 



10% 



(54) 



The above expression is an upper limit on the magnitude of Af^f k given any choice of k. For a mi which is large enough 
to be non- negligible, we need only require an accuracy in its measurement of Ami = 0.1 and a measurement accuracy for 
B GGG and C GG of 10% in order to have flf k < O(10~ 4 ), which is safely negligible by a factor of 10 compared to the minimum 
measurement error e™™ of the lensing bispectrum. As discussed by Zhang, this level of accuracy can likely be accomplished by 
direct measurement of rrii under the approximation C\l' « C GG (see Eq. 7) if the lensing contamination C// < 10%. However, 
if the II contamination is greater than 10% of the lensing signal, more detailed methods must be employed to achieve a great 
enough accuracy in the ra» measurement for it to be safely negligible, some of which are discussed by Zhang (2010a). 



4.4 Non-Gaussianity and galaxy bias 

We are only interested in the bispectrum due to the non-linear evolution of gravitational clustering and the associated 
intrinsic alignment contamination, leaving the accurate estimation of the bispectrum due to primordial non-Gaussianity to 
other works. Equation 13 should then include a term B^ G ° which must be separately accounted for. Similarly, the relation 
between the 3D matter bispectrum and 3D galaxy bispectrum depends on non-Gaussianity beyond the scale dependent 
correction b 1 (z) — > bi(z) + Ab(k, z) used in relating the 3D matter power spectrum to the 3D galaxy power spectrum (Jeong 
& Komatsu 2009). Equation 22 must also include the contributions by non-Gaussianity in the term B^ G ° and from the 
trispectrum which we have previously neglected. 

The full expression including all non-Gaussian contributions is given in Appendix B of Jeong & Komatsu (2009). However, 
from Figs. 10-14 of Jeong & Komatsu (2009), it is clear that if we avoid very stretched or elongated triangle shapes that are 
very sensitive to non-Gaussianity, at the scales of interest in a lensing survey (10 2 < i < 10 4 ), the total contribution to the 
relation by non-Guassianity as a fraction of the non-linear term is less than 10% for Jnl = 40 and qnl = 10 4 . If we accept 
the smaller values of /nl = 4 and qnl = 100, this fractional contribution is less than 1%. Thus we can safely ignore the 
contribution of the non-Gaussianity as a source of error to the relation since it is expected to be on the order of the minimum 
GG measurement error and less than the systematic error discussed in Sec. 4.2. Future work will better constrain and model 
the effects of non-Gaussianity, thus allowing its effect to be fully accounted for in the GGI self-calibration. 

The linear galaxy bias is discussed by Zhang (2010a), and the error induced by the expected uncertainty in its measure- 
ment in the GI self-calibration is demonstrated to be negligible compared to other sources of error. The linear and non-linear 
galaxy bias terms can be measured simultaneously by using the approach of Fry (1994). Using measurements of C" and 
B 933 , we extract the bias information from the relationship 



B 333 (£ 1 ,£2,£ 3 ) « (b\) 3 BZr m + (b\)% (C™ m (£i)C™ m (£ 2 ) + C™ m (£ 2 )C™ m (£ 3 ) + C™ m (h)^™ (£ 3 )) , (55) 

where C mm and B mmm are the angular matter power and bispectrum, weighted identically to galaxies. Equation 55 is 
the analog to Eq. 22, which includes an intrinsic alignment component. The matter power and bispectrum can be tightly 
constrained by CMB measurements and then evolved, given a cosmology, to low redshift to predict C mm and B rnmm . 

Both measurement error in B 333 and uncertainties in the predictions of C mm and B mmm will affect the measurement 
of the linear and non-linear galaxy bias parameters. The non-linear galaxy bias is more difficult to constrain precisely than 
the linear galaxy bias, with typical measured and expected uncertainties in its measurement of up to A\b\/ {b l 2 ) 2 ] « 0.5 at 
la confidence (LSST Science Collaborations and LSST Project 2009; Simpson et al. 2011). We use the estimate of Zhang 
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Figure 5. The fractional error ), = AB IGG / B IGG in Eq. 27 is shown for equilateral triangles (£ = £\ = £2 = £3) due to uncertainties 



ijk 

in the linear and non-linear galaxy bias. e^ k is typically less than 2% except for large scales. The measurement error this induces in the 
final measurement of B GGG , Af.' 6 ? — 

' J ilk 



*-3 



\fijk' ls re< iuced by the factor f^ k ^ 1. It is thus typically negligible when compared to the 



minimum measurement error e 



ijk 



111 



B GGG ( Fig 3) 



(2010a) for the measurement error in the linear galaxy bias 



to plot in Fig. 5 the fractional error e\X = AB IGG / B IGG in Eq. 27 of even a large uncertainty for the non- linear galaxy bias 
of Ab 2 ~ 0.5. We find that eS is generally less than 2% except for large scales. The measurement error this induces in the 



Afej 

,IGG 



eAef, 



s k I/ 



C33 



(56) 



e ijkfijk- Even for very large /y fe = 1, Af- jk is typically comparable to or less 

jGGG 



final measurement of B GGG is then A/^ 

than the minimum measurement error e™j k n in B GGG (Fig. 3). For a typical f'j^, we would expect it to be entirely negligible 
for all bin choices. 

The only real limitation which comes from the galaxy bias is then the scale to which it can be applied in the non-linear 
regime. Recent work (Simpson et al. 2011) has shown that the scale down to which the bias model we have employed is 
accurate can be extended to k = 0.5, which corresponds to I m 1000 at the median redshift of an LSST-like survey. Future 
work may extend this range further, but for now this places an approximate upper limit on the £ at which the self-calibration 
can function to a high degree of accuracy. In the future, a more robust bias model could be chosen for the very highly non- 
linear regime to extend this limit with relative ease, as it will alter only the form of Eq. 27 and the resulting performance 
calculations, while the method of extracting B\f® remains unchanged. 



4.5 Other sources of uncertainty 

The GGI self-calibration requires the calculation of Wijk and Qs, which include the cosmology-dependent lensing kernel. 
This introduces an uncertainty due to the measurement of fl m and the distance-redshift relation. However, we expect this 
uncertainty to be negligible when compared to other dominant sources of error in the GGI self-calibration. Komatsu et al. 
(2011) have measured Q m to 5% accuracy, and new measurements are expected to constrain Q m to 1-2%. The distance-redshift 
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relation will also be constrained to 1% by baryon acoustic oscillations and supernovae (Albrecht et al. 2006). We expect that 
given these constraints, any uncertainty introduced by the lensing kernel will only affect the GGI self-calibration at the percent 
level, which is negligible compared to the expected systematic error Sfijk of Eq. 51. An iterative approach can also be applied, 
where a set of initial cosmological parameters is chosen as above and used for the 2- and 3-point self-calibration, from which 
new (improved) parameter constraints can be calculated and applied again until the interactive process converges. 

Similarly, we have used an approximate fitting formula derived from perturbation theory by Scoccimarro & Couchman 
(2001) for the bispectrum in our error estimations. This is only expected to be accurate to within 15% when compared to 
N-body simulations for the lensing bispectrum. We thus expect uncertainty due to the calculation of the bispectrum to be 
dominant when compared to errors associated with the power spectrum calculation. A more accurate approach to modelling the 
bispectrum and the effects of intrinsic alignment will provide more accurate estimates of the GGI self-calibration performance, 
which we leave to a later work. 

Catastrophic photo-z error also affects the GGI self-calibration through the assumed galaxy distribution. We assume 
a Gaussian photo-z PDF in our numerical calculations, but observed photo-z PDFs generally have non-negligible outliers. 
This affects the GGI self-calibration through the calculation of Qz and the relationship between B IGG and B 1 " . However, 
these effects are suppressed due to both numerator and denominator being affected in similar ways. The effect can be further 
decreased by better photo-z PDF template estimates and better calibration of photo-z errors, and we expect the GGI self- 
calibration to ultimately be safe from non- negligible degradation due to catastrophic photo-z errors. 

The relationship between B IGG and B 1 " depends upon our assumption of a deterministic galaxy bias, which is not 
perfectly accurate in real galaxy distributions. This could cause both random and systematic error in the GGI self-calibration. 
A true quantification of this effect is beyond the scope of this paper, as the possible correlation between stochasticity and 
intrinsic alignment is not well understood. However, Baldauf et al. (2010) has shown that it is possible to suppress the galaxy 
stochasticity to the 1% level in some cases, which allows that the effect of stochasticity in the GGI self-calibration could 
ultimately be limited to the percent level, which would be safely negligible compared to other sources of error. 

4.6 Summary of residual errors 

There are three regimes under which the performance of the GGI self-calibration can be summarised. These are defined by 
the magnitude of the GGI contamination as represented by f-j k . The first is where the ggl correlation is too small to detect 
in B^ 2 \ with f(j k ^ fij* k esh - If the intrinsic alignment cannot be detected in B^\ the GGI self-calibration is not applicable. 
This generally means that the GGI contamination is also negligible when compared to e*™jL n , the minimum statistical error in 
the lensing bispectrum, and there is no need to correct for it. 

If fijk > ftj'k esh S then the GGI contamination to the lensing bispectrum is likely not negligible, and it must be corrected 
for. The GGI self-calibration is now able to detect and calculate the GGI correlation. In the second regime, where Af^ k esh > 
€ ijkfijky trie statistical error induced by measurement error in the estimator B'? 9 is dominant. As shown in Fig. 3, this 

error is generally negligible when compared to e™™, and so in this regime, the GGI self-calibration should perform at the 
statistical limit of the lensing survey. 

Finally, where Af^ k reah < e*?£flj k , the systematic error Sf ijk = <^ k flj k due to the relationship between B\ GG and B\?? 
in Eq. 27 is dominant. In the case where t s ^ k < t'ij k 1 '/fij k , e^jk 1 is still dominant. Otherwise the GGI self-calibration can 
suppress the GGI contamination by a factor of 5-10 or more for all but a few adjacent redshift bin choices. In this case, other 
complementary techniques could be employed to further reduce the GGI contamination down to the statistical limit for the 
lensing survey. 

In the 2-point correlations, one such case has been explored by Zhang (2010b), but such studies of the 3-point intrinsic 
alignment are left to be done. Zhang, Pen & Bernstein (2010) combines the GI self-calibration with a photo-z self-calibration 
to better protect the GI self-calibration against catastrophic photo-z effects. Both methods are possible because the GI and 
GGI self-calibration uses primarily those correlations in one redshift bin to estimate the intrinsic alignment, while Zhang 
(2010b); Zhang, Pen & Bernstein (2010) use those correlations between redshift bins. As first mentioned in Sec. 1, others 
have also used information between redshift bins to calibrate the intrinsic alignment contamination in the 2- and 3-point 
correlations (Okumura T., Jing 2009; Kirk, Bridle & Schneider 2010; Joachimi & Schneider 2008, 2009; Shi, Joachimi & 
Schneider 2010; Joachimi & Bridle 2010). Such techniques for the 3-point intrinsic alignment correlations should eventually 
complement the GGI self-calibration for improved reductions in the contamination by the intrinsic alignment in the cosmic 
shear signal, but much work is left to be done. 



5 CONCLUSION 

The GGG bispectrum has been shown to be strongly contaminated by the 3-point intrinsic alignment correlations. While the 
III and Gil correlations can be neglected by considering only the cross-correlation bispectrum between three different redshift 
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bins, the GGI correlation remains a contaminant. Zhang (2010a) first proposed the self-calibration technique in order to 
calculate and remove the 2-point GI contamination from the GG power spectrum. In this work we verify the performance of 
the GI self-calibration technique, and expand the self-calibration to the 3-point correlations, proposing the GGI self-calibration 
technique to calculate and remove the GGI correlation from the GGG bispectrum. 

We first establish the estimator B*? 9 to extract the ggl correlation from the galaxy ellipticity-density-density measurement 
for a photo-z galaxy sample. We show that this estimator is expected to be generally applicable to weak lensing surveys and 
reduces to the simple extraction method for spectroscopic galaxy samples at low photo-z error. We then develop a relation 
between the GGI and ggl bispectra using the linear and non-linear galaxy bias to relate the galaxy density and cosmic shear 
measurements. This allows us to calculate and remove the GGI correlation from the GGG bispectrum. While this method 
is in principle applicable to all i and triangle shapes, we do note some modest restrictions in section 4.4 on very elongated 
triangles due to the effects of non-Gaussianity and at very non-linear scales due to limitations in the understanding of the 
galaxy bias model used. 

We quantify the performance of the GGI self-calibration technique for a typical weak-lensing survey, using anticipated 
parameters for the LSST as an example case. The residual statistical error due to measurement uncertainty in the estimator 
Biii i s sh° wn to be generally negligible when compared to the minimum measurement error in the lensing bispectrum. By 
considering the systematic error introduced by the relationship between B\f h a and B{f?, we show that for galaxy triplets 
with bins which are not adjacent, \tijk\ < 0.1. For these bin choices, the intrinsic alignment contamination can be suppressed 
by a factor of 10 or greater. In most cases where two or three bins are adjacent, \tijk\ < 0.2, which allows for a suppression 
in the contamination by a factor of 5. In only a few of the cases where all three bins are adjacent is \tijk\ > 0.2, and even 
in these cases we expect a suppression in the contamination by a factor of 3 or more. This will potentially allow the GGI 
self-calibration to reduce the GGI correlation to the statistical limit of the lensing survey, as discussed in Sec. 4.6. 

These results are insensitive to the original intrinsic alignment contamination, such that for any fijk esh < fijk < 1> the 
GGI self-calibration will reduce the GGI contamination down to survey limits or by a factor of 5-10 or greater, whichever is 
less, for all but a few adjacent redshift bin triplets. This is only slightly reduced from the GI self-calibration, where for any 
fthresh < < QI self-calibration reduces the GI contamination down to survey limits or by a factor of 10 or greater, 

whichever is less. We thus expect the GGI self-calibration to perform near the level of the GI self-calibration, and together 
they promise to be an efficient technique to isolate both the 2- and 3-point intrinsic alignment signals from the cosmic shear 
signal. 
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APPENDIX A: CALCULATION OF COEFFICIENTS IN AB, 



IGG 

m 



Upon evaluating the sum and taking the Fourier transform of Eq. 43, each of the products of the correlations have a numerical 
coefficient due to the restrictions on redshift ordering. Many, however, are identical due to symmetries. The calculation of 
the unique coefficients a-h in Eqs. 45 & 46 are summarised here. The first coefficient is trivial, due to products with no 
noise correlations or correlations like (8 a n v ), which are themselves orientation dependent. We then calculate for a term like 

{8 a Sx){8p8 ll ){K n K v ) 

AT-6 

(1 _ Q a Y2 E I> 5 "^ - Q3)(3Sa m , - Qs) « 1. (Al) 



.2 

ck/37 A/if 



For terms like {5 a S x ){Sp8 tM ){K- f K v ) oc 8 a \, which include one galaxy density noise correlation 

a = (l NP Q 3) 2 E E( 3S ^ - Qs)(3S aM , - Qs) « 1 + gnVp - (A2) 
For terms like (S a Sx)(SpSfj,){K^ oc 5^, which include one convergence noise correlation 

b = (I-Q3) 2 ^ E( 3 ^ - Q3)(35 AM7 - Q 3 ) » 1 + 5(1 _ 4 Qa)2 . (A3) 

a/37 
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For terms like {S a 8 X ){Sp 8^ }(k 7 k„) oc SaxSp^, which include two galaxy density noise correlations 



— (i _ Q 3 )2 E E( 3S «^ - Q3)(3^ 7 - Q 3 ) - 1 - 4(1 _ 1 Q3)2 - (A4) 

For terms like (8™ $\}{5i38fj,}(Ky k^} oc 8 a \8~,v, which include one galaxy density noise correlation and one convergence noise 
correlation 



d = (IZ%)2 E E( 35 ^ ^ 0a)(35. M7 - Qs) « 1 + 4(1 _ Qa)2 . (A5) 
For terms like (5^ Sx)(Sp 8^}{k^ k„) oc 6 a \6p,j.5 lv , which include only noise correlations 



For terms like (<5 a <5,\)(K-y<S/i)(«i/<5/9), which include two correlations with the orientation dependence described in Eq. 44 



A7 — ° 9 

(A6) 



5Q§ - Qa(3 + 13Q 2 + 2g_|) + 5Q 2 (1 + 2Q 2 ) , , : ; 



80Q 2 (1-Q 2 ) 2 (1-Q 3 ) 2 

For terms like (8™ 8 x ){K^8^}{K v 8fs) oc 5 a A, which include two correlations with the orientation dependence described in Eq. 
44 and one galaxy density noise correlation 



E E( 3 ^ - Qs)(3S a »„ Qs) 1 - ( + §f ) \ ( (i#^y + ff 



3 + 5Qg - Q 3 (3 + 13Q 2 + 2Qj) + 15Q| 

80Q|(1-Q 2 ) 2 (1-Q 3 ) 2 ' 1 ' 

Finally, for terms like (8„ 8 x )(K, 1 8 tl )(I v 8f i ) oc 8 a \, which include one correlation with the orientation dependence described in 
Eq. 44 and one galaxy density noise correlation 

h = ( i %.A2 E E( 3S <^ - Q3)(3S QM „ - Qs)\ (t^t + |f 

20Ql - 5Q 3 (5 + 6Q 2 ) + 6(1 + 6Q 2 ) . . 

80Q 2 (1-Q 2 )(1-Q 3 ) 2 ' 1 j 

For typical values Q 2 = 1/2 and Q 3 = 2/5, these coefficients are 

. . 

d = « 4.5 

36 

59 

e = y «6.6 
/ = 1 

. - s-» 

h = y « 1.5. (A10) 



© 0000 RAS, MNRAS 000, 1-22 



GGI self- calibration 21 



REFERENCES 

Acquaviva V., Hajian A., Spergel D., Das S., 2008, PRD, 78, 043514 

Albrecht A., et al., 2006, Report of the Dark Energy Task Force, arXiv:astro-ph/0609591 

Bacon D.J., Refregier A.R., Ellis R.S., 2000, MNRAS, 318, 625 

Bacon D.J., Refregier A., Clowe D., Ellis R.S., 2001, MNRAS, 325, 1065 

Baldauf T., Smith R., Seljak U., Mandelbaum R., 2010, PRD, 81, 063531 

Bean R., Tangmatitham M., 2010, PRD, 81, 083534 

Bernstein CM., 2009, ApJ, 695, 652 

Bernstein CM., Jarvis M., 2002, AJ, 123, 583 

Blazek J., McQuinn M., Seljak U., 2011, JCAP, 05, 010 

Bridle S., King L., 2007, New J. Phys. 9, 444 

Brown M.L., Taylor A.N., Hambly N.C, Dye S., 2002, MNRAS, 333, 501 

Brown M.L., Taylor A.N., Bacon D.J., Gray M.E., Dye S., Meisenheimer K., Wolf C, 2003, MNRAS, 341, 100 
Capozzicllo S., Cardone V.F., Troisi A., 2006, PRD, 73, 104019 
Catelan P., Kamionkowski M., Blandford R.D., 2001, MNRAS, 320, L7 
Cooray A., Sheth R., 2002, Phys. Rept. 372, 1 

Crittenden R.C, Natarajan P., Pen U.-L., Theuns T., 2001, ApJ, 559, 552 
Croft R., Metzler C, 2000, ApJ, 545, 561 

Daniel S., Caldwell R., Cooray A., Melchiorri A., 2008, PRD, 77, 103513 

Daniel S., Linder E., Smith T., Caldwell R., Corray A., Leauthaud A., Lombriser L., 2010, PRD, 80, 123508 
Dossett J., Moldenhauer J., Ishak M., 2011, PRD, 84, 023012 
Eisenstein D.J., Hu W., Tegmark M., 1999, ApJ, 518, 2 

Erben T., Van Waerbeke L., Bertin E., Mellier Y., Schneider P., 2001, A&A, 366, 717 

Faltenbacher A., Li C, White S.D.M., Jing Y.P., Mao S., Wang J., 2009, Res. Astron. & Astrophys., 9, 41 

Fry J., 1994, PRL, 73, 215 

Fry J.N., Gaztanaga E., 1993, ApJ, 413, 447 

Fu X., Wu P., Yu H., 2009, PLB, 677, 12 

Heavens A., Refregier A., Heymans C, 2000, MNRAS, 319, 649 
Heymans C, Heavens A., 2003, MNRAS, 339, 711 

Heymans C, Brown M., Heavens A., Meisenheimer K., Taylor A., Wolf C, 2004, MNRAS, 347, 895 
Heymans C, White M., Heavens A., Vale C, Van Waerbeke L., 2006, MNRAS, 371, 750 
Hirata CM., Seljak U., 2003, MNRAS, 343, 459 
Hirata CM., Seljak U., 2003, PRD, 67, 43001 
Hirata CM., Seljak U., 2004, PRD, 70, 063526 

Hirata CM., Mandelbaum R., Ishak M., Seljak U., Nichol R., Pimbblet K.A., Ross N.P., Wake D., 2007, MNRAS, 381, 1197 

Hoekstra H., Yee H.K.C., Gladders M.D., Barrientos L.F., Hall P.B., Infante L., 2002, ApJ, 72, 55 

Hu W., 2002, PRD, 65, 023003 

Hu W., Tegmark M., 1999, ApJL, 514, L65 

Huterer D., Linder E., 2007, PRD, 75, 023519 

Ishak M., Dossett J., 2009, PRD, 80, 043004 

Ishak M., Hirata CM., McDonald P., Seljak U., 2004, PRD, 69, 083514 
Ishak M., Upadhye A., Spergel D., 2006, PRD, 74, 043513 

Jarvis M., Bernstein CM., Fischer P., Smith D., Jain B., Tyson J. A., Wittman D., 2003, AJ, 125, 1014 

Jing Y.P., 2002, MNRAS, 335, 89 

Jeong D., Komatsu E., 2009, ApJ, 703, 1230 

Joachimi B., Bridle S., 2010, A&A, 523, Al 

Joachimi B., Schneider P., 2008, A&A, 488, 829 

Joachimi B., Schneider P., 2009, A&A, 507, 105 

Joachimi B., Schneider P., 2010, A&A, 517, A4 

Joachimi B., Mandelbaum R., Abdalla F., Bridle S., 2010, A&A, 527, A26 

Joudaki S., Cooray A., Holz D.E., 2009, PRD, 80, 023003, 

King L., 2005, A&A, 441, 47 

King L., Schneider P., 2002, A&A, 396, 411 

King L., Schneider P., 2003, A&A, 398, 23 

Kirk D., Bridle S., Schneider M., 2010, MNRAS, 408, 1502 

Komatsu E., et al., 2011, ApJS, 192, 18 

© 0000 RAS, MNRAS 000, 1-22 



22 Troxel & Ishak 



Krause E., Hirata CM., 2011, MNRAS, 410, 2730 
Linder E., Calm R., 2007, Astropart. Phys. 28, 481 

LSST Science Collaborations and LSST Project, 2009, LSST Science Book, Version 2.0, arXiv:astro-ph/0912.0201, 

http: / /www. lsst . org/lsst / scibook / 
Mandelbaum R., Hirata CM., Ishak M., Seljak U., Brinkmann J., 2006, MNRAS, 367, 611 
Massey R., Refregier A., Bacon D., Ellis R., 2005, MNRAS, 359, 1277 
Okumura T., Jing Y.P., 2009, ApJL, 694, L83 

Pen U.-L., Lu T., Van Waerbeke L., Mellier Y., 2003, MNRAS, 346, 994 

Refregier A., 2003, Ann. Rev. A&A, 41, 645 

Rhodes J., Refregier A., Groth E.J., 2001, ApJL, 552, L85 

Schmidt F., 2008, PRD, 78, 043002 

Schrabback T. et al., 2010, A&A, 516, A63 

Scoccimarro R., Couchman H., 2001, MNRAS, 325, 1312 

Semboloni E., Heymans C, Van Waerbeke L., Schneider P., 2008, MNRAS, 388, 991 

Semboloni E., Schrabback T., Van Waerbeke L., Vafaei S., Hartlap J., Hilbert S., 2010, MNRAS, 410, 143 

Simpson F., James J., Heavens A., Heymans C, 2011, arXiv:astro-ph/1107.5169 

Shi X., Joachimi B., Schneider P., 2010, A&A, 523, A60 

Song Y.S., 2005, PRD, 71, 024026 

Takada M., Jain B., 2003, MNRAS, 340, 580 

Takada B., Jain M., 2004, MNRAS, 348, 897 

Takada M., White M., 2004, ApJL, 601, LI 

Thomas S., Abdalla F., Weller J., 2009, MNRAS, 395, 197 

Toreno I., Semboloni E., Schrabback T., 2010, A&A, 530, A68 

Vafaei S., Lu T., Van Waerbeke L., Semboloni E., Heymans C, Pen U.-L., 2010, Astropart. Phys. 32, 340 
Van Waerbeke L., Mellier Y., 2003, ArXiv Astrophysics e-prints, astro-ph/0305089 
Van Waerbeke L. et al., 2000, A&A, 358, 30 

Van Waerbeke L., Mellier Y., Pell R., Pen U.-L., McCracken H.J., Jain B., 2002, A&A, 393, 369 
Zaldarriaga M., Spergel D.N., Seljak U., 1997, ApJ, 488, 1 
Zhang P., 2010, ApJ, 720, 1090 
Zhang P., 2010, MNRAS, 406, L95 

Zhang P.J., Liguori M., Bean R., Dodelson S., 2007, PRL, 99, 141302 
Zhang P., Pen U.-L., Bernstein C, 2010, MNRAS, 405, 359 
Zhao H., Bacon D.J., Taylor A.N., Home K., 2006, MNRAS, 368, 171 
Zhao G., Pogosian L., Silvestri A., Zylberberg J., 2009, PRD, 79, 083513 



© 0000 RAS, MNRAS 000, 1-22 



